Effect of mutations on the p53 IRES RNA structure Implications for de-regulation of the synthesis of p53 isoforms

Earlier we have demonstrated the presence of internal ribosome entry site (IRES) within tumor suppressor p53 mRNA. Here we have mapped the putative secondary structure of p53-IRES RNA using information from chemical probing and nuclease mapping experiments. Additionally, the secondary structure of the IRES element of the wild-type RNA was compared with cancer-derived silent mutant p53 RNAs. These mutations might result in the conformational alterations of p53-IRES RNAs. The results also indicate decreased IRES activities of the mutants as compared to wild-type RNA. Further, it was observed that some of the cytoplasmic trans-acting factors, critical for enhancing IRES function, were unable to bind mutant RNAs as efficiently as to wild-type. Our results suggest that hnRNP C1/C2 binds to p53-IRES and siRNA mediated partial silencing of hnRNP C1/C2 showed appreciable decrease in IRES function and consequent decrease in the level of the corresponding p53 isoform. Interestingly mutant p53 IRES showed lesser binding with hnRNP C1/C2 protein. Finally, upon doxorubicin treatment, the mutant RNAs were unable to show enhanced p53 synthesis to similar extent compared to wild type. Taken together, these observations suggest that mutations occurring in the p53 IRES might have profound implications for de-regulation of its expression and activity.

Cardioviral RNA structure logo analysis: entropy, correlations, and prediction

In recent years, there has been an increased number of sequenced RNAs leading to the development of new RNA databases. Thus, predicting RNA structure from multiple alignments is an important issue to understand its function. Since RNA secondary structures are often conserved in evolution, developing methods to identify covariate sites in an alignment can be essential for discovering structural elements. Structure Logo is a technique established on the basis of entropy and mutual information measured to analyze RNA sequences from an alignment. We proposed an efficient Structure Logo approach to analyze conservations and correlations in a set of Cardioviral RNA sequences. The entropy and mutual information content were measured to examine the conservations and correlations, respectively. The conserved secondary structure motifs were predicted on the basis of the conservation and correlation analyses. Our predictive motifs were similar to the ones observed in the viral RNA structure database, and the correlations between bases also corresponded to the secondary structure in the database.

Computational Methods for RNA Structure Validation and Improvement.

With increasing recognition of the roles RNA molecules and RNA/protein complexes play in an unexpected variety of biological processes, understanding of RNA structure-function relationships is of high current importance. To make clean biological interpretations from three-dimensional structures, it is imperative to have high-quality, accurate RNA crystal structures available, and the community has thoroughly embraced that goal. However, due to the many degrees of freedom inherent in RNA structure (especially for the backbone), it is a significant challenge to succeed in building accurate experimental models for RNA structures. This chapter describes the tools and techniques our research group and our collaborators have developed over the years to help RNA structural biologists both evaluate and achieve better accuracy. Expert analysis of large, high-resolution, quality-conscious RNA datasets provides the fundamental information that enables automated methods for robust and efficient error diagnosis in validating RNA structures at all resolutions. The even more crucial goal of correcting the diagnosed outliers has steadily developed toward highly effective, computationally based techniques. Automation enables solving complex issues in large RNA structures, but cannot circumvent the need for thoughtful examination of local details, and so we also provide some guidance for interpreting and acting on the results of current structure validation for RNA.

Mutational interference mapping experiment (MIME) for studying RNA structure and function

RNA regulates many biological processes; however, identifying functional RNA sequences and structures is complex and time-consuming. We introduce a method, mutational interference mapping experiment (MIME), to identify, at single-nucleotide resolution, the primary sequence and secondary structures of an RNA molecule that are crucial for its function. MIME is based on random mutagenesis of the RNA target followed by functional selection and next-generation sequencing. Our analytical approach allows the recovery of quantitative binding parameters and permits the identification of base-pairing partners directly from the sequencing data. We used this method to map the binding site of the human immunodeficiency virus-1 (HIV-1) Pr55(Gag) protein on the viral genomic RNA in vitro, and showed that, by analyzing permitted base-pairing patterns, we could model RNA structure motifs that are crucial for protein binding.

Altered 3′-terminal RNA structure in phage Qβ adapted to host factor-less Escherichia coli

The RNA phage Qβ requires for the replication of its genome an RNA binding protein called Qβ host factor or Hfq protein. Our previous results suggested that this protein mediates the access of replicase to the 3′-end of the Qβ plus strand RNA. Here we report the results of an evolutionary experiment in which phage Qβ was adapted to an Escherichia coli Q13 host strain with an inactivated host factor (hfq) gene. This strain initially produced phage at a titer ≈10,000-fold lower than the wild-type strain and with minute plaque morphology, but after 12 growth cycles, phage titer and plaque size had evolved to levels near those of the wild-type host. RNAs isolated from adapted Qβ mutants were efficient templates for replicase without host factor in vitro. Electron microscopy showed that mutant RNAs, in contrast to wild-type RNA, efficiently interacted with replicase at the 3′-end in the absence of host factor. The same set of four mutations in the 3′-terminal third of the genome was found in several independently evolved phage clones. One mutation disrupts the base pairing of the 3′-terminal CCCoh sequence, suggesting that the host factor stimulates activity of the wild-type RNA template by melting out its 3′-end.

An RNA structure involved in feedback regulation of splicing and of translation is critical for biological fitness.

While studies of the regulation of gene expression have generally concerned qualitative changes in the selection or the level of expression of a gene, much of the regulation that occurs within a cell involves the continuous subtle optimization of the levels of proteins used in macromolecular complexes. An example is the biosynthesis of the ribosome, in which equimolar amounts of nearly 80 ribosomal proteins must be supplied by the cytoplasm to the nucleolus. We have found that the transcript of one of the ribosomal protein genes of Saccharomyces cerevisiae, RPL32, participates in such fine tuning. Sequences from exon I of the RPL32 transcript interact with nucleotides from the intron to form a structure that binds L32 to regulate splicing. In the spliced transcript, the same sequences interact with nucleotides from exon II to form a structure that binds L32 to regulate translation, thus providing two levels of autoregulation. We now show, by using a sensitive cocultivation assay, that these RNA structures and their interaction with L32 play a role in the fitness of the cell. The change of a single nucleotide within the 5' leader of the RPL32 transcript, which abolishes the site for L32 binding, leads to detectably slower growth and to eventual loss of the mutant strain from the culture. Experiments designed to assess independently the regulation of splicing and the regulation of translation are presented. These observations demonstrate that, in evolutionary terms, subtle regulatory compensations can be critical. The change in structure of an RNA, due to alteration of just one noncoding nucleotide, can spell the difference between biological success and failure.

A satellite RNA of barley yellow dwarf virus contains a novel hammerhead structure in the self-cleavage domain

An RNA molecule with properties of a satellite RNA was found in an isolate of barley yellow dwarf virus (BYDV), RPV serotype. It is 322 nucleotides long, single-stranded, and does not hybridize to the viral genome. Dimers of the RNA, which presumably represent replicative intermediates, were able to self-cleave into monomers. In vitro transcripts from cDNA clones were capable of self-cleavage in both the plus (encapsidated) and minus orientations. The sequence flanking the minus strand cleavage site contained a consensus " hammerhead" structure, similar to those found in other self-cleaving satellite RNAs. Although related to the hammerhead structure, sequences flanking the plus strand termini showed differences from the consensus and may be folded into a different structure containing a pseudoknot. © 1991.

Translation of the flavivirus Kunjin NS3 gene in cis but not its RNA sequence or secondary structure is essential for efficient RNA packaging

Our previous studies using trans-complementation analysis of Kunjin virus (KUN) full-length cDNA clones harboring in-frame deletions in the NS3 gene demonstrated the inability of these defective complemented RNAs to be packaged into virus particles (W. J. Liu, P. L. Sedlak, N. Kondratieva, and A. A. Khromykh, J. Virol. 76:10766-10775). In this study we aimed to establish whether this requirement for NS3 in RNA packaging is determined by the secondary RNA structure of the NS3 gene or by the essential role of the translated NS3 gene product. Multiple silent mutations of three computer-predicted stable RNA structures in the NS3 coding region of KUN replicon RNA aimed at disrupting RNA secondary structure without affecting amino acid sequence did not affect RNA replication and packaging into virus-like particles in the packaging cell line, thus demonstrating that the predicted conserved RNA structures in the NS3 gene do not play a role in RNA replication and/or packaging. In contrast, double frameshift mutations in the NS3 coding region of full-length KUN RNA, producing scrambled NS3 protein but retaining secondary RNA structure, resulted in the loss of ability of these defective RNAs to be packaged into virus particles in complementation experiments in KUN replicon-expressing cells. Furthermore, the more robust complementation-packaging system based on established stable cell lines producing large amounts of complemented replicating NS3-deficient replicon RNAs and infection with KUN virus to provide structural proteins also failed to detect any secreted virus-like particles containing packaged NS3-deficient replicon RNAs. These results have now firmly established the requirement of KUN NS3 protein translated in cis for genome packaging into virus particles.

miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star

On the role of RNA silencing in the pathogenicity and evolution of viroids and viral satellites

Viroids and most viral satellites have small, noncoding, and highly structured RNA genomes. How they cause disease symptoms without encoding proteins and why they have characteristic secondary structures are two longstanding questions. Recent studies have shown that both viroids and satellites are capable of inducing RNA silencing, suggesting a possible role of this mechanism in the pathology and evolution of these subviral RNAs. Here we show that preventing RNA silencing in tobacco, using a silencing suppressor, greatly reduces the symptoms caused by the Y satellite of cucumber mosaic virus. Furthermore, tomato plants expressing hairpin RNA, derived from potato spindle tuber viroid, developed symptoms similar to those of potato spindle tuber viroid infection. These results provide evidence suggesting that viroids and satellites cause disease symptoms by directing RNA silencing against physiologically important host genes. We also show that viroid and satellite RNAs are significantly resistant to RNA silencing-mediated degradation, suggesting that RNA silencing is an important selection pressure shaping the evolution of the secondary structures of these pathogens.

The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes

In Arabidopsis thaliana (Arabidopsis), DICER-LIKE1 (DCL1) functions together with the double-stranded RNA binding protein (dsRBP), DRB1, to process microRNAs (miRNAs) from their precursor transcripts prior to their transfer to the RNA-induced silencing complex (RISC). miRNA-loaded RISC directs RNA silencing of cognate mRNAs via ARGONAUTE1 (AGO1)-catalyzed cleavage. Short interefering RNAs (siRNAs) are processed from viral-derived or transgene-encoded molecules of doublestranded RNA (dsRNA) by the DCL/dsRBP partnership, DCL4/DRB4, and are also loaded to AGO1-catalyzed RISC for cleavage of complementary mRNAs. Here, we use an artificial miRNA (amiRNA) technology, transiently expressed in Nicotiana benthamiana, to produce a series of amiRNA duplexes with differing intermolecular thermostabilities at the 5′ end of duplex strands. Analyses of amiRNA duplex strand accumulation and target transcript expression revealed that strand selection (amiRNA and amiRNA*) is directed by asymmetric thermostability of the duplex termini. The duplex strand possessing a lower 59 thermostability was preferentially retained by RISC to guide mRNA cleavage of the corresponding target transgene. In addition, analysis of endogenous miRNA duplex strand accumulation in Arabidopsis drb1 and drb2345 mutant plants revealed that DRB1 dictates strand selection, presumably by directional loading of the miRNA duplex onto RISC for passenger strand degradation. Bioinformatic and Northern blot analyses of DCL4/DRB4-dependent small RNAs (miRNAs and siRNAs) revealed that small RNAs produced by this DCL/dsRBP combination do not conform to the same terminal thermostability rules as those governing DCL1/DRB1-processed miRNAs. This suggests that small RNA processing in the DCL1/DRB1-directed miRNA and DCL4/DRB4-directed sRNA biogenesis pathways operates via different mechanisms.

Small RNA sequencing of potato leafroll virus-infected plants reveals an additional subgenomic RNA encoding a sequence-specific RNA-binding protein

Potato leafroll virus (PLRV) is a positive-strand RNA virus that generates subgenomic RNAs (sgRNA) for expression of 3' proximal genes. Small RNA (sRNA) sequencing and mapping of the PLRV-derived sRNAs revealed coverage of the entire viral genome with the exception of four distinctive gaps. Remarkably, these gaps mapped to areas of PLRV genome with extensive secondary structures, such as the internal ribosome entry site and 5' transcriptional start site of sgRNA1 and sgRNA2. The last gap mapped to ~500. nt from the 3' terminus of PLRV genome and suggested the possible presence of an additional sgRNA for PLRV. Quantitative real-time PCR and northern blot analysis confirmed the expression of sgRNA3 and subsequent analyses placed its 5' transcriptional start site at position 5347 of PLRV genome. A regulatory role is proposed for the PLRV sgRNA3 as it encodes for an RNA-binding protein with specificity to the 5' of PLRV genomic RNA. © 2013.

RNA silencing platforms in plants

Since the discovery of RNAi, its mechanism in plants and animals has been intensively studied, widely exploited as a research tool, and used for a number of potential commercial applications. In this article, we discuss the platforms for delivering RNAi in plants. We provide a brief background to these platforms and concentrate on discussing the more recent advances, comparing the RNAi technologies used in plants with those used in animals, and trying to predict the ways in which RNAi technologies may further develop. © 2005 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.